Heat transport and diffusion in a canonical model of a relativistic gas 
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Relativistic transport phenomena are important from both theoretical and practical point of view. 
Accordingly, hydrodynamics of relativistic gas has been extensively studied theoretically. Here, we 
introduce a three-dimensional canonical model of hard-sphere relativistic gas which allows us to 
impose appropriate temperature gradient along a given direction maintaining the system in a non- 
equilibrium steady state. We use such a numerical laboratory to study the appropriateness of the 
so-called first order (Chapman-Enskog) relativistic hydrodynamics by calculating various transport 
coefficients. Our numerical results are consistent with predictions of such a theory for a wide range of 
temperatures. Our results are somewhat surprising since such linear theories are not consistent with 
the fundamental assumption of the special theory of relativity {v < c). We therefore seek to explain 
such results by studying the appropriateness of diffusive transport in the relativistic gas, comparing 
our results with that of a classical gas. We find that the relativistic correction (constraint) in the 
hydrodynamic limit amounts to small negligible corrections, thus indicating the validity of the linear 
approximation in near equilibrium transport phenomena. 

PACS numbers: 47.75.-l-f,51.10.-|-y,02.70.Ns,95.30.Lz 



I. INTRODUCTION 

Owing to its diverse applications in the hydrodynamic 
description of astrophysical phenomena [l|-|j], heavy ion 
coUision experiments and hot plasmas generated in col- 
Uders relativistic kinetic theory has attracted con- 

siderable attention ever since it was formulated. The ex- 
tension of the theory to non-equilibrium situations also 
covers different problems ranging from formation of the 
universe to engineering applications such as pulsed laser 
cutting and welding or high speed machining processes 
[13) [13] ■ Perhaps the central issue in this field is the 
relation between thermodynamic forces and fluxes in a 
system out of equilibrium. From a phenomenological 
point of view, in hydrodynamic regime and not too far 
from equilibrium, fluxes are linearly related to thermody- 
namic forces with the so called transport coefficients as 
proportionality factors. The theories that assume such 
linear relations are known as first- order theories and are 
based on the hypothesis of local equilibrium in the hy- 
drodynamic regime. A familiar example is the Fourier 
law of heat transport that relates heat current to the 
temperature gradient. The objection raised against such 
linear laws is that they lead to diffusion equations imply- 
ing an unbounded propagation of fluctuations in sharp 
contrast to the first assumption of the special theory of 
relativity. Various suggestions, including Israel-Stewart 
second-order formulation |12Hl5j and several more ex- 
tensions p^ - fTsj are proposed in the literature to remove 
the non-causality and instabilities associated with first 
order theories. The price to be paid, however, is the 
introduction of complicated higher-order terms into the 
formulations as also noted in [T^, [2^ . It is not our goal 
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here to add a new proposal or even judge the existing 
theories. Here, we rather propose to study how well (or 
badly) do first-order theories capture the general features 
of relativistic transport phenomena. 

Our approach is motivated by the need to find the 
simplest approximate model given both the significance 
of such problems and the complicated existing alterna- 
tives, mentioned above. We therefore propose and study 
a canonical three-dimensional (3d) model of a relativis- 
tic hard-sphere gas in contact with pre-defined tempera- 
tures. We are therefore able to impose temperature gra- 
dients on our system and use such a realistic numerical 
laboratory to calculate various transport coefficients in 
the intermediate and highly relativistic regimes. We then 
compare our results with the theoretical values obtained 
from Chapman-Enskog (CE) approximation to relativis- 
tic transport equations which constitutes a well-known 
example of a linear first order theory. We find good 
agreement between our simulation results and theoretical 
predictions. We then seek to explain such agreements by 
studying diffusive transport in classical versus relativistic 
dynamics. In Section |lTl the main features of relativistic 
kinetic theory as well as CE approximation are reviewed. 
Section IIIII describes the details of our numerical model 
in order to simulate a relativistic hard-sphere gas in con- 
tact with a heat bath. In Section ITVl the results of our 
simulations and a comparison with previous analytical 
calculations are presented. Finally, we close with con- 
cluding remarks in Section IVl 



II. RELATIVISTIC KINETIC THEORY 

In a statistical description of many-body systems, the 
one-particle distribution function is a key ingredient. It 
can be interpreted as giving the average number of par- 
ticles with a certain momentum at each space-time point 
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and its explicit expression is obtaine d by solving a kinetic 
equation called transport equation [2l|, [22] ■ This equa- 
tion which was first derived by Boltzmann for a classical 
gaseous system, gives the rate of change of the distri- 
bution function in time and space due to the particle 
interactions. Generally speaking, the form of transport 
equation is based upon three assumptions: (i) only binary 
collision are considered, (ii) the spatiotemporal changes 
of distribution function on microscopic scales are negli- 
gible and (iii) the "assumption of molecular chaos" (or 
"Stosszahlansatz" [2^) holds. The last one is a statistical 
assumption which says no correlations exists between the 
colliding particles; that is, the average number of binary 
collisions is proportional to the product of their distribu- 
tion function. It is therefore evident how the microscopic 
structures of the system are reflected in the distribution 
function and why it can relate microscopic features to 
equilibrium and non-equilibrium macroscopic properties. 

The relativistic counterpart of the classical transport 
equation is obtained following the same reasoning, with 
same assumptions as above, but within the framework of 
relativistic theory 24 26]. 



p^dJ{x,p)^C[f{x,p)] 



(1) 



in which f{x,p) is the distribution function, x and p are 
respectively position and momentum four- vectors, and 



d^, = d/dx^ 



'^dt,V). The collision function, which 



includes the particle interactions has the form 



C[f] 



1 



d^pi d^p' dVi 



2 J p{ p - p^ 
[f.f[W{p',p'^\p,p,) ~ //ilF(p,Pi|p',p'i)], (2) 



where W denotes the transition rate and is determined 
by the underlying dynamics [23 |. The abbreviations / 
and /i are used for f{x,p) and f{x,pi) and the primes 
refer to quantities after collision. The solution of the 
homogeneous transport equation, i.e. C[f] = 0, leads to 
the equilibrium distribution. 



/ 



exp[-/3^p^] 
(27r;i)3 ■' 



(3) 



in which (27rfi,) denotes Planck constant and (3^ = 
-j^j^U^ is the temperature four- vector. The hydrody- 
namic velocity four-vector, = ^{u){c,u) with 7(u) = 
^/l — {u/cY as Lorentz factor, is defined so that UfJJ^ = 
(? . Eq. ([3]) which is best known as Jiittner distribution, is 
the relativistic generalization of the Maxwell-Boltzmann 
distribution and its validity has been a source of long 
lasting debates in the physics community [27l - l32| . Hav- 
ing this distribution in hand, the particle four-flow in 
equilibrium is given by 



(4) 



and the equilibrium particle density is defined as 
n := c~^J{^Uf^ \24^ . Similarly, the equilibrium energy- 



momentum tensor can be defined as 

(^/^PVexphM, (5) 



and the energy density as ne := c~'^T^'^U ^Ujj. 

When the system is not too far from equilibrium, the 
transport equation may be solved assuming the distribu- 
tion function is a function of local hydrodynamic vari- 
ables, e.g. particle density, energy density and hydro- 
dynamic velocity, as well as their gradients. The sys- 
tematic method to obtain such solutions was worked out 
independently by Chapman ^33] and Enskog (s^]. Their 
approach was based on an expansion of distribution func- 
tion around local equilibrium distribution, /(°\ with the 
ratio of mean free path and a typical macroscopic length 
as an expansion parameter, here denoted by -q: 



(1) 



i2f(2) 



(6) 



The main interest of CE method lies in the first correc- 
tion term, which leads to linear laws for the transport 
phenomena. Of course, this approximation can only be 
considered accurate in the hydrodynamic regime where 
the expansion parameter is small. 

The classical CE approximation can be used in rela- 
tivistic kinetic theory without major modifications [sH, 
[36| . In the first approximation the one-particle distribu- 
tion is written as: 



f{x,p)^f(°Hx,p)[l + <f{x,p)] 



(7) 



where f^°'>{x,p) is the local relativistic distribution func- 
tion. Following the standard kinetic theory, substitution 
of the appropriate form of $, as a function of hydrody- 
namic variables and their gradients, into the relativistic 
heat flux expression would giv e the relativistic general- 
ization of Fourier's law f2J, [33, 



T 



(8) 



in which, P is pressure, e and h = e 
and enthalpy per particle, and A'^'^ := 



Pn , are energy 



'WW' 



IS 

the projection operator with = diag(l, — 1, — 1, — 1) 
the metric tensor. It can be shown that the heat trans- 
port coefficient to first order, in CE approximation for a 
three dimensional hard-sphere gas is [23, [13] 



3 cfcB z^j^|(z)[7/(7-l)]^ 
327r a (z^ + 2)K2{2z) + 5zK3{2z) 



(9) 



where z = and Kn's are the Bessel functions of the 
second kind [38J, 7 = cp/cy, and a is the differential 
cross section. Using the equation of state of an ideal 
relativistic gas, P — nk^T, to convert pressure to particle 
density and defining as A^'^d^, Eq.® can be written 
in the familiar form 



TT- 



(10) 
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where one defines 

Ltt = AT(1-^), (11a) 

Ltu = ^-f^ (lib) 

to respectively distinguish the contributions of energy 
flux and particle flux in the total heat flow [STj . 

III. DEFINITION OF THE MODEL 

In (3lj | we introduced a two dimensional model of a 
relativistic hard-disk gas which can be used as a numer- 
ical laboratory to investigate various prop erties of a rel- 
ativistic gas in the low density regime [sil . |32| . [soj . Here, 
we use a three dimensional generalization of the previ- 
ous model in which the particles are spheres rather than 
disks. The interaction cross-section associated with hard- 
disk (or sphere in 3d) model is independent of the energy 
and of the scattering angle. This makes it a suitable 
choice to simulate high energy hadrons since experimen- 
tal evidences show that the total cross sections for mas- 
sive hadron scattering are more or less constant in the 
energy range of interest in relativistic kinetic theory . 
For instance, the total cross-section for nucleon-nucleon 
scattering is roughly constant from an energy of 2GeV 
on [i^l- This is the typical energy for the hadron era of 
the big bang theory |24 |. 

The system we study here, consists of TV particles of 
equal rest mass m which are constrained to move in a 
box of volume V = L^LyLz- The spherical particles 
move in straight lines at constant speed and change their 
momenta instantaneously when they touch at distance 5. 
The binary collision of particles in three dimension in- 
volves six degrees of freedom (three degrees for each par- 
ticle) which can be solved deterministically if we limit 
ourselves to head-to-head collisions. In such collisions, 
the force is exerted along the line connecting the centers, 
'I'ij — fi ^ Tj, and thus the two components of parti- 
cles' momenta perpendicular to Tij remain unchanged. 
The parallel components, however, change according to 
energy- momentum conservation laws as below (29l , l3l| : 

P'i,\\ = liVcmfl'iVcn.ei - (1 + Wcrn)Pi,||]> 

e- = l{vcrnf[{l + vlJe, - 2v,^p^\l (12) 

where hatted quantities refer to momenta after collision 
and Vcm — {p i,\\~^Pj.\\) / {Ei-\-Ej) is the collision invariant, 
relativistic center-of-mass velocity of the two particles. 
With the same rules for particle j, a deterministic, time- 
reversible canonical transformation at each collision is 
defined. 

A hard-sphere model such as above avoids the general 
problems associated with interacting relativistic particle 
systems [4l|. Thus it can be successfully applied to sim- 
ulate a relativistic gas to shed light on challenging ques- 
tions like how thermostatistical properties change when 



the gas is viewed by a moving observer as compared 
to a rest observer [3l| or what happens if the measure- 
ments are performed with respect to laboratory time [sij , 
proper time [32] or on the light-cone ji^]. These issues 
were studied in microcanonical ensemble where any type 
of interaction with environment is absent. In this article 
we introduce thermal interactions and/or temperature 
gradients to the system to perform simulations in the 
canonical ensemble or to study out-of-equilibrium prop- 
erties. 

The heat bath we have used here, is a stochastic 
thcrmo-bath [s^ |4^, 13] which interacts locally, i.e. 
through the walls, with the particles inside the box. Ef- 
fectively, it is a random number generator that simulates 
the exchange of energy between an ideal relativistic gas 
and the canonical relativistic system at a predefined tem- 
perature T. Whenever a particle impinges on the bound- 
aries, the components of its velocity are updated. More 
specifically, we remove the incident particle and introduce 
a new particle with a new random energy value e' with 
the following rates, while keeping its position unchanged: 

W{e{v) e'{v')) = 1 for Vy' > 

W{e{v) ^ e'{v')) = ^ for V,, < V, (13) 

This choice of rates is consistent with the condition of 
detailed balance 45]. In order to get the correct veloc- 
ity distribution, the probability is defined, according to 
Jiittner distribution, as 

oc 7(w)''+2 exp[-;3m7(w)], (14) 

in which v represents the magnitude of the velocity, d is 
the space dimension {d — i here), and the temperature 
of the bath is introduced through the parameter P = 
l/k^T [sil, H^. Such thermalization process allows the 
particle's energy to be exchanged if the new energy is 
more probable (Vv' > Vv) while it gives the system a 
chance to deviate from the desired distribution {Vv' < 
Vv)- Accepting energies with reduced probability allows 
the system to sample fluctuations away from the lowest 
energy configuration and thus gives the correct Jiittner 
distribution. 

In our numerical simulations, the walls normal to the 
x-direction are diathermal and other walls are periodic. 
Therefore, the sign of the j:-component of the updated 
velocity of the particle is manually chosen to make sure 
that the incident particle reflects back into the box. Also, 
the components of the updated velocity should be chosen 
appropriately so that its magnitude respects the main 
assumption of relativity, \v\ < c. 

Numerical results indicate that the effective tempera- 
ture obtained using the above stochastic bath does not 
exactly meet the nominal value of bath's temperature, 
especially when the walls are kept in different tempera- 
tures. This observation which is also reported in other 
non- relativistic transport models [HjH^] would not affect 
our general results; i.e., the existence of a well-defined 
stationary state, linear temperature profile and the va- 
lidity of local thermal equilibrium [39[ . 
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FIG. 1. (color online). Constant temperature profile T{x) 
for a canonical simulation. TV = 256 particles were put in a 
box of size LxLyL^ with its diathermal walls kept at the same 
temperature. The system is divided into equally-spaced verti- 
cal bins and the temperature for each bin has been calculated 
using Eq. (|17|) . Inset shows the measured velocity distribution 
of particles which is well fitted to the Jiittner function 
(solid red line) with the expected temperature. 

IV. RESULTS 

In the following sections the results of our simulations 
in canonical ensemble and out-of-equilibrium states are 
presented. In all simulations, the values of rest mass m, 
particles' diameter S, Boltzmann constant fee, and light 
velocity c, are all set to unity. To satisfy the condition 
of zero net momentum in the rest frame, half the par- 
ticles are given random initial velocities and the second 
half are chosen to have the same velocities with opposite 
signs as the first half. The total momentum remains zero 
(on the average) after the system is allowed to interact 
with the heat baths and thus the assumption of zero hy- 
drodynamic velocity is valid. To obtain the velocity dis- 
tributions, we let the system equilibrate (typically after 
lO^iV collisions) and simultaneously measure velocities 
of all particles with respect to a rest frame at a given 
instant of time t. To collect more data we may either re- 
peat the above procedure for different initial conditions 
or perform measurements at several well-separated time 
instances. 



A. Canonical ensemble 

The model described above lead to a well-defined equi- 
librium state in the canonical ensemble. To show this we 
have measured the velocity distribution of the particles 
with respect to a rest observer and laboratory time t [i^ . 
The result shown in Fig |T] clearly indicates that the par- 
ticles have reached a state consistent with the prediction 
of equilibrium relativistic statistics, i.e., the Jiittner ve- 
locity distribution and a constant temperature profile. 
To obtain the temperature profile we typically divide our 
system into M equal bins of volume L^LyLz/M . Each 
bin is large enough to allow computation of (local) rele- 




FIG. 2. (color online), (a) Shows changes in temperature 
of leftmost and rightmost bins immediately after the system 
is exposed to a temperature gradient. In the resulting non- 
equilibrium stationary state, the velocity distribution of the 
middle bin of the system is well fitted to the Jiittner func- 
tion with the expected temperature. The linear behavior of 
temperature profile (c) and the resulting gradient in particle 
density (d) are also shown. 

vant thermodynamic variable. For instance, the average 
number of particles in each bin is computed by 

f I ft N 

"fe=/ y,Kx-x,{t'))dt'dx, (15) 

J Ax * Jo 

in which, Xi(t') is the position of i-th particle at time t' 
and Ax is the region that contains position of the /c*'* 
bin, i.e., [xk — ^,Xk + The same technique is used 
to obtain the average energy, e^, in each bin: 

f 1 /■* ^ 

/ T y,Six-x,it'))e,it')dt'dx, (16) 

J Ax t Jo -^^ 

The temperature profile, T^, is then computed as a func- 
tion of rife and using the following equation [sil li^ : 

— = — +m— — — , (17) 

rik Pk K2{mj3k) 

which is written for d = 2>. 



B. Out-of-equilibrium studies 

We now turn to the out of equilibrium situation in 
which the vertical left and right walls are kept at differ- 
ent temperatures. Results of simulation with N = 256 
particles and box of size = 45, Ly = Lz = 10 indicate 
that after typically lO^A^ collisions, a non-equilibrium 
stationary state is obtained. In this state the average 
value of thermodynamic properties remain constant and 
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FIG. 3. (color online). The numerical (dots) and theoreti- 
cal (lines) heat transport coefficients of a relativistic gas as a 
function of temperature. The coefficients Ltt (upper branch) 
and Ltu (lower branch), are normalized to the low tempera- 
ture limit Ao — 75fcB kBT / rmv / {I28a) . The nominal values 
of temperature difference in all simulations is AT — 0.5. Note 
that in the relativistic high temperature regime the contribu- 
tion of particle flux is comparable to the contribution of en- 
ergy flow as predicted by the relativistic theory of irreversible 
processes to linear order. 



the condition of local thermal equilibrium is satisfied (see 
Fig|5|). That is, for small temperature gradients, in every 
macroscopically infinitesimal volume element, thermody- 
namic variables can be defined in the usual way and are 
related to each other through the relations which hold for 
the system at equilibrium. To show this, we have divided 
the system into nine bins. The temperature profile, Tk, is 
computed as discussed in the previous section. Figmja) 
shows how the temperature of different bins (here the 
leftmost and rightmost bins) change after the introduc- 
tion of different baths until a non-equilibrium stationary 
state is achieved. In this state a linear temperature pro- 
file is obtained (Figl^c)) which indicates that VT{x) is 
constant throughout the system. The temperature gra- 
dient imposed, on the other hand, causes a gradient in 
particle density (Figl^Jd)) as expected from Onsagcr re- 
lations Additionally, the velocity distribution for a 
typical bin (see FigHJ^b)) is well fitted to Jiittner distribu- 
tion with the corresponding temperature from FiglSfc). 
These conditions indicate that the obtained stationary 
state is locally equilibrated and thus can be used to study 
transport properties of relativistic ideal gases. 



C. Heat transport coefficients 

Our main concern in this section is to study the heat 
transport properties of the hard-sphere relativistic gas. 
We therefore propose to verify the validity of linear ap- 
proximation (i.e. Eq pU)) ). by directly calculating the 
heat flux in the presence of a temperature gradient and 
thus calculating the relevant transport coefficients and 
comparing them with theoretical prediction Eq. pT|) . In 
order to numerically obtain the value of A in this equa- 
tion, we impose a fixed temperature gradient, as de- 
scribed above, and directly measure the energy fiow in 
our simulations using the following formula: 

Je ^ ^ (j^^-'^^y (18) 

where (. . .) denotes time averaging. Note that the system 
is simulated in the rest frame, i.e the average total mo- 
mentum of the system in all directions and thus the av- 
erage hydrodynamic velocity is zero. Therefore, the ori- 
gin of energy flow is purely chaotic, that is = Jq [Slj . 
Since the temperature gradient is imposed in the x direc- 
tion the only nonzero component of heat current is Jg . 
Other numerically measured quantities are VT, Vn, T, 
and h. As stated in Sec. IIIIl the temperature of walls do 
not (exactly) coincides with the nominal temperature of 
the heat baths. Therefore, the effective temperature gra- 
dient obtained in simulations, VT*, as well as numerical 
value of Vn is substituted in Eqs. (|TOl) . The values of 
temperature, T, and enthalpy, h, are mean values evalu- 
ated for the middle bin. Due to linear behavior of temper- 
ature profile, the temperature of the middle bin can also 
be expressed as the average of the corresponding values 
in the rightmost and leftmost bins. It is now straightfor- 
ward to calculate heat transport coefficient. A, and conse- 
quently the coefficients Ltt and Ltu, defined in Ea.(|lip. 
Figure [3] shows the results obtained by averaging over 10^ 
realizations. The coefficients are normaliz ed to the low 
temperature limit of Eq. dH), Aq = ^-^\f^ H- As 
indicated, the behavior of the numerically obtained trans- 
port coefficients agree well with theoretical predictions of 
relativistic CE approximation, (Eqs.® and (ITT|) ). Our 
results therefore show the validity of such linear approx- 
imation for a wide range of temperatures well into the 
relativistic regime. 

It is worth emphasizing that the dynamics of the hard- 
sphere model we have used here fully respects the as- 
sumptions of relativity, in the sense that the velocity of 
particles as well as propagation speed of fluctuations in 
thermodynamic properties such as energy, heat, density, 
etc. are all less than the maximum allowed velocity, c. 
Under such circumstances, the agreement between nu- 
merical results and theoretical predictions of relativistic 
CE approximation indicates that this linear theory is ca- 
pable of obtaining heat transport coefficients of a rela- 
tivistic gas. The validity of CE approximation, however. 
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is limited by the ratio of mean free path to the character- 
istic length scale (or relaxation time scale to the dynam- 
ical time scale) of the system which forms the expansion 
parameter 77 in Eq.®. This parameter should be small 
enough to let us ignore nonlinear terms. It is therefore 
expected that CE fails when superfluid helium, rarefied 
gasses, or nanoscale problems are studied. This classi- 
cally accepted criterion also holds in relativistic regime. 
However, it is not straightforward to numerically inves- 
tigate the validity of higher order terms via the present 
(low density) model. In this model we have neglected is- 
sues such as the possibility of radius contraction of mov- 
ing spheres, the possibility of multiple (non-binary) colli- 
sions, and the importance of particle interactions beyond 
hard sphere potential. We are not aware of a determin- 
istic MD model concerning all theses issues though in 
some kinetic models, Boltzmann transport equation has 
been solved for a dense ultra-relativistic gas considering 
binary (and including) three-body interactions [7.: ^] • 

Another important result is the contribution of parti- 
cle flow in the flux of heat (see FiglH). As is well-known, 
heat flow is due to the diffusive transport of energy in 
a non-relativistic limit and the contribution of particle 
flow, Ltu, in the flux of heat is well negligible. How- 
ever, as can be seen in FigEl the contribution of particle 
flow becomes an increasingly significant part of heat fiux 
as one moves towards the relativistic, high temperature 
regime. This purely relativistic effect add more subtlety 
to the theory and thus could lead to interesting nontrivial 
consequences when other thermodynamic are considered 
(See below for more discussions.). 



D. Diffusive transport of heat 



In Section IIV CI we showed that the heat flux in a 
relativistic system is well described by the linear ap- 
proximation Eq. (jlO[) . Since no net flux of particles 
is present in the non-equilibrium stationary state (i.e. 
Jji = LnT^j^ — Lnn^-^ = 0), it is possible to write 
Eq. ([TU]) in the closed form 



(19) 



in which k = [L'^-^x — LTTLnn)/ (TLnti) is the thermal 
conductivity. Combining this linear law with the con- 
servation of energy-momentum tensor {df/T^^ — 0) and 
using the relation T"*'^ = cyT, which is valid for both 
low and high temperature limits, results in the following 
diffusion equation for propagation of heat: 



dT{r,t) 



(20) 



(a) 



0.008 
^ 0.006 

0.004 
a, 0.002 




t=5, classic 
t=5, relativistic 
t=20, classic 
t=20, relativistic 



-20 




X 



20 




Relativistic 
Classic 



(b) 



200 
t 



FIG. 4. (color online). Part (a) shows short time and long 
time behavior of relativistic propagators obtained for a 3d 
gaseous system with particles of diameter 5 = 1, number den- 
sity n — 0.064. Both systems are initiated with same ini- 
tial conditions. The corresponding mean free time and tem- 
perature of classical and relativistic gas are computed to be 
MFTci. = 3.57, MFTrei. = 3.55, Td. = 0.33 and Trei. = 2.66. 
Part (b) shows the diffusion coefficient as a function of time 
that fits to the curve D{t) = 1.85 — for classical data and 
agrees well with the analytically predicted value = 1.89. 
Note that the relativistic coeflicient is always less than the 
classical value due to the constraints of relativity. 



is obtained for density: 

dn{r, t) 
dt 



(21) 



The diffusive propagation of heat / density is known to 
be inconsistent with the basic assumption of relativity as 
it allows superluminal transport. The underlying mech- 
anism for such macroscopic relaxation process is the ran- 
dom walk nature of Brownian particles also described by 
a diffusion equation whose solution is given by propaga- 
tor in the d-dimensional space. 



p{r,t\ro,0) 



1 



exp(- 



(22) 



where Pt — ^i^/cvi T^'^ denotes energy-momentum ten- 
sor and cv is the molar heat capacity. On the other hand, 
if the system is put between two particle reservoirs with 
different chemical potentials, the same type of equation 



where T) is the appropriate diffusion constant, and Tq 
is the position of the Brownian particle at t = 0. The 
propagator give the probability of finding the particle at 
time t in position r. The nonzero probability at arbitrary 
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long distances for short times is the basis of the objec- 
tion to such diffusive properties within relativistic theo- 
ries. While such an objection is clearly justified, the ques- 
tion we ask is that under typical relativistic parameters 
how differently does a Brownian particle move as com- 
pared to its classical counterpart. We therefore calculate 
such propagators directly by monitoring a "tagged" par- 
ticle. In order to simulate identical situations, the initial 
conditions (i.e. initial positions and velocities) are cho- 
sen to be the same in both classical and relativistic sys- 
tems. Fig.Ul^a) shows the propagators in the x-direction. 
As expected, the classical propagator fits to a Gaussian 
which is basically different from relativistic propagator 
in the sense that the latter is restricted to the region 
a; — a;o| < c(< — to) while the former allows an unbounded 
propagation of particles (or heat). This fundamental dif- 
ference is the origin of objections against diffusive trans- 
port of particles (or heat) and consequently against the 
so called first order theories. Given such discrepancy, the 
question is why the relativistic first order theory of CE is 
good enough when diffusion or heat transport coefficients 
of a relativistic system are considered? The reason traces 
back to the hydrodynamic limit in which the long time 
(compared to mean free time) behavior of the system is 
considered. 

As Fig. m^a) indicates, in this limit the difference of 
classical and relativistic propagators diminishes so that 
the classical Gaussian propagator (namely, the solution 
of diffusion equation) becomes an almost exact approxi- 
mation of the relativistic propagator. Consequently, the 
resulting relativistic diffusion coefficient shows no ma- 
jor distinction from its classical value as indicated in 
Fig. HKb). In this figure the classical as well as rel- 
ativistic diffusion coefhcient are plotted as a function 
of time. The transient behavior seen between ballistic 
(< (Ar)^ >cx t"^) and diffusive (< (Ar)^ >cx t) regimes 
is a sign of correlations ( or memory effects) in the sys- 
tem which is well studied in the literature |54l | . It can be 
shown that the diffusion coefficient of a classical hard- 
sphere gas fits to the curve D{t) = V — which, af- 
ter sufficiently long time, saturates to the constant co- 
efficient predicted for a simple uncorrelated model (i.e., 
T> = 3 /(8nS'^) ^/ k-p,T / rmr with n denoting number den- 
sity [5^.). Also, note that the value of relativistic diffu- 
sion is always less than the classical one due to the well 
known relativistic constraints. The important point to 
be noted here is the similarity of classical and relativistic 
coefficients which clearly supports the appropriateness of 
first order theories and the consequent diffusion equation 
in the hydrodynamic limit. 

Here, we would like to make an important point. De- 



spite the fact that we have argued that diffusion equation 
is a fairly good approximation in the hydrodynamic limit, 
there is still a fundamental difference between the clas- 
sical and relativistic propagators in Fig Surely, they 
look very similar, but nevertheless one is a Gaussian and 
the other is not. The non-Gaussian nature of the rel- 
ativistic propagator is due to the relativistic constraint 
\v\ < c which effectively introduces some sort of corre- 
lation/memory in the random walk process thus making 
it a non-markovian process [i^, [s^. Such correlations 
may have important consequences as far as cooperative 
phenomena are considered. If such phenomena exist we 
expect them to become more pronounced in the high tem- 
perature (relativistic) regime as such correlations become 
more pronounced. This is potentially very interesting 
since cooperative phenomena in non-interacting systems 
are rare and thus unexpected results. A closer study of 
this issue is out of the scope of this manuscript and de- 
serves a separate publication [56i] . 

V. CONCLUDING REMARKS 

In this article we have presented results of numerical 
simulation of a three-dimensional canonical generaliza- 
tion of our previous model of a relativistic gas. We have 
used such a numerical laboratory in order to study trans- 
port properties in the presence of a temperature gradi- 
ent. We use parameter regimes where local equilibrium is 
shown to be well-established. Our main result is that var- 
ious transport coefficients are consistent with the CE (lin- 
ear) approximation to the relativistic transport equation. 
We observe such consistency for a wide range of tempera- 
ture well into the high relativistic regime. Our results are 
somewhat surprising because of the well-known inconsis- 
tency between first order theories and the fundamental 
assumption of special relativity. We have therefore pro- 
posed to look into transport propagators both for the 
classical as well as highly relativistic gas in order to find 
an explanation for the observed accuracy of the linear 
approximation. We find that in the long time hydro- 
dynamic limit, the kinetic description of such diffusive 
processes are almost identical, with practically negligi- 
ble difference between the classical and the relativistic 
regime. The Gaussian nature of such propagators and 
the resulting diffusion coefficient is much the same as the 
relativistic one despite the fact that the latter respects 
the constraints of special relativity. Our results provide 
a basis for justification of applicability of first order rel- 
ativistic transport theories despite their well-known lim- 
itations. 
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